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Abstract 



An exactly solvable position-dependent mass Schrodinger equation in two dimensions, de- 
picting a particle moving in a semi-infinite layer, is re-examined in the light of recent 
theories describing superintegrable two-dimensional systems with integrals of motion that 
are quadratic functions of the momenta. To get the energy spectrum a quadratic algebra 
approach is used together with a realization in terms of deformed parafermionic oscillator 
operators. In this process, the importance of supplementing algebraic considerations with 
a proper treatment of boundary conditions for selecting physical wavefunctions is stressed. 
Some new results for matrix elements are derived. Finally, the two-dimensional model is 
extended to two integrable and exactly solvable (but not superintegrable) models in three 
dimensions, depicting a particle in a semi-infinite parallelepipedal or cylindrical channel, 
respectively. 
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1 Introduction 



Quantum mechanical systems with a position-dependent (effective) mass (PDM) have at- 
tracted a lot of attention and inspired intense research activites during recent years. They 
are indeed very useful in the study of many physical problems, such as electronic properties 
of semiconductors pQ and quantum dots [2], nuclei [3], quantum liquids [I], 3 He clusters [5], 
metal clusters [Bj, etc. 

Looking for exact solutions of the Schrodinger equation with a PDM has become an 
interesting research topic because such solutions may provide a conceptual understanding 
of some physical phenomena, as well as a testing ground for some approximation schemes. 
Although mostly one-dimensional equations have been considered up to now, several works 
have recently paid attention to <i-dimensional problems [TJ, [HI El El El 12] • 

In [9] (henceforth referred to as I and whose equations will be quoted by their number 
preceded by I), we have analyzed the problem of d- dimensional PDM Schrodinger equations 
in the framework of first-order intertwining operators and shown that with a pair (H, Hi) 
of intertwined Hamiltonians we can associate another pair (R, Ri) of second-order partial 
differential operators related to the same intertwining operator and such that H (resp. 
Hi) commutes with R (resp. Ri). In the context of supersymmetric quantum mechanics 
(SUSYQM) based on an sl(l/l) superalgebra, R and Ri can be interpreted as SUSY part- 
ners, while H and Hi are related to the Casimir operator of a larger gl(l/l) superalgebra. 

In the same work, we have also applied our general theory to an explicit example, 
depicting a particle moving in a two-dimensional semi-infinite layer. This model may be of 
interest in the study of quantum wires with an abrupt termination in an environment that 
can be modelled by a dependence of the carrier effective mass on the position. It illustrates 
the influence of a uniformity breaking in a quantum channel on the production of bound 
states, as it was previously observed in the case of a quantum dot or a bend p2] . 

From a theoretical viewpoint, our model has proved interesting too because it is solvable 
in two different ways: by separation of variables in the corresponding Schrodinger equation 
or employing SUSYQM and shape-invariance techniques [13] . The former method relies 
upon the existence of an integral of motion L, while, as above-mentioned, the latter is 
based on the use of R. In other words, the three second-order partial differential operators 
H, L and R form a set of algebraically independent integrals of motion, which means that 
the system is superintegrable. 
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Let us recall that in classical mechanics |15j . an integrable system on a <i-dimensional 
manifold is a system which has d functionally independent (globally defined) integrals of 
motion in involution (including the Hamiltonian) . Any system with more that d functionally 
independent integrals of motion is called superintegrable. It is maximally superintegrable if 
it admits the maximum number 2d— 1 of integrals of motion. The latter form a complete set 
so that any other integral of motion can be expressed in terms of them. In particular, the 
Poisson bracket of any two basic integrals, being again a constant of motion, can be written 
as a (in general) nonlinear function of them. Such results can be extended to quantum 
mechanics [16] , so that for quantum counterparts of maximally superintegrable systems we 
get (in general) nonlinear associative algebras of algebraically independent observables, all 
of them commuting with H. 

The simplest case corresponds to the class of two-dimensional superintegrable systems 
with integrals of motion that are linear and quadratic functions of the momenta. The 
study and classification of such systems, dating back to the 19th century and revived in the 
1960ties [IT], have recently been the subject of intense research activites and substantial 
progress has been made in this area (see [TH HH [20l [2H [22l [23l [2H [251 (2SI [271 [2B] and 
references quoted therein). In particular, it has been shown that their integrals of motion 
generate a quadratic Poisson algebra (in the classical case) or a quadratic associative algebra 
(in the quantum one) with a Casimir of sixth degree in the momenta and the general 
form of these algebras has been uncovered [221 I2H]- Algebras of this kind have many 
similarities to the quadratic Racah algebra QR(3) (a special case of the quadratic Askey- 
Wilson algebra QAW(3)) [19]. They actually coincide with QR(3) whenever one of their 
parameters vanishes. The eigenvalues and eigenf unctions of the superintegrable system 
Hamiltonian can be found from the finite-dimensional irreducible representations of these 
algebras. The latter can be determined by a ladder-operator method [191 [2U] or through a 
realization [2H [22] in terms of (generalized) deformed parafermionic operators [29] , which 
finite-dimensional version of deformed oscillator operators [3U] . 

Since our two-dimensional PDM model belongs to this class of superintegrable systems, 
it is interesting to analyze it in the light of such topical and innovative theories. This is one 
of the purposes of the present paper, which will therefore provide us with a third method for 
solving the PDM Schrodinger equation. In such a process, we will insist on the necessity of 
supplementing algebraic calculations with a proper treatment of the wavefunction boundary 
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conditions imposed by the physics of the problem — a point that is not always highlighted 
enough. 

The other purpose of the present paper is to free ourselves from the restriction to a two- 
dimensional model. We actually plan to show that an abrupt termination of a quantum 
channel can also be mimicked by some three-dimensional exactly solvable models. 

This paper is organized as follows. In Section 2, the two-dimensional PDM model of I 
is briefly reviewed and some important comments on its mathematical structure are made 
in conjunction with the physics of the problem. The quadratic algebra approach to such 
a model is then detailed in Section 3. Two three-dimensional extensions of the model are 
presented in Section 4. Finally, Section 5 contains the conclusion. 

2 Exactly solvable and superintegrable PDM model 
in a two-dimensional semi-infinite layer 

In I we considered a particle moving in a two-dimensional semi-infinite layer of width ix/q, 
parallel to the x-axis and with impenetrable barriers at the boundaries. The variables x, y 
vary in the domain 

7T 7T 

D : < x < oo, < y < — , (2.1) 

2q y 2q' V ' 

and the wavefunctions have to satisfy the conditions 



7T 



#),j/) = 0, ^k±- =0. (2.2) 

The mass of the particle is m{x) = moM(x), where the dimensionless function M{x) is 
given by 

M(x) = sech 2 qx. (2.3) 
In units wherein h = 2itlq = 1, the Hamiltonian of the model can be written as 

where 

V£jf(x) = —q 2 cosh 2 qx + q 2 k(k — 1) csch 2 qx (2.5) 

is an effective potential including terms depending on the ambiguity parameters (see 
Eq. (12.3)). In f)2.5p . the constant k is assumed positive and we have set an irrelevant 
additive constant v to zero. 
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Both the operators 
and 



(2.6) 



— cosh 2 qx sin 2 qyd^ + 2 sinh qx cosh qx sin qy cos qy d 2 xy — sinh 2 ga; cos 2 qy <9 2 

+ q sinh qx cosh gx(l — 4 sin 2 gy)^ + g(l + 4 sinh 2 ga;) sin qy cos gy9 y 

+ g 2 (sinh 2 qx — sin 2 qy — 3 sinh 2 gx sin 2 gy) — g 2 /c(l + csch 2 qx sin 2 gy) 

+ g 2 /c 2 csch 2 qx sin 2 gy, (2.7) 



where 



^0)t _ _ gQg^ q X s j n ^ _|_ q X cos qyQ y _q gj^ ^ s j n ^ 

— g/c csch qx sin gy, 

77^ = cosh ga; sin gy d x — sinh ga; cos qy d y + q sinh qx sin gy 

— g/c cschga; singy, 



(2.8) 
(2.9) 



commute with H^ k \ although not with one another. Hence one may diagonalize either 
and L or and simultaneously. This leads to two alternative bases for the 
Hamiltonian eigenfunctions, corresponding to the eigenvalues 



EP = q\N + 2){N + 2k + 1), N = 0, 1, 2, . . . , 



with degeneracies 



deg(TV) 



N' 
~2 



+ 1, 



(2.10) 



(2.11) 



where [iV/2] stands for the integer part of N/2. 

The first basis is obtained by separating the variables x, y in the PDM Schrodinger 
equation and its members, associated with the eigenvalues (/ + l) 2 g 2 of L, read 



i> { n}(x,y) = €i(x)xi(y), n,/ = o,i,2,..., 



(2.12) 



with N = 2n + Z, 



=A^y(tanhga:) fe (sechga;) i+2 Pr 2 '' +lj 



(1 - 2tanh 2 qx), 



(2.13) 



/f cos[(Z + l) ? y] for/ = 0,2,4,..., 

xi(y)={ V ( 2 - 14 ) 

sin[(Z + l)gy] for/ = 1,3,5,..., 

and TV^y given in Eq. (13.18). 

The second basis, resulting from the intertwining relation 

v (k) H (k) = H (k) v (k) } Jjik) = H (k+1) + 2q 2^ (2 15) 

and its Hermitian conjugate, can be built by successive applications of operators of type 

(x, y) = ^ ^(»+i)t . . . ,(*+-Dt*g^) y) , (2.16) 

on functions W^^J o (x, y), annihilated by r/( fe+l/ ) and given in Eqs. (13.28), (13.32) and (13.34). 
In (12.161) . iV runs over 0, 2, 4,. . . , N or iV — 1, according to whether N is even or odd, 
while v, defined by v = N — N , determines the R( k > eigenvalue 

r ( fc ) = q 2 v(v + 2k), v = 0,1,2,.... (2.17) 

Although an explicit expression of the normalization coefficient J\f^ No is easily obtained 
(see Eq. (13.41)), this is not the case for ^f^ No (x,y) (except for some low values of iV and 
No), nor for the expansion of *nn ( x i v) ^ ne ^ rs ^ basis eigenfunctions ip^J(x, y), which 
is given by rather awkward formulas (see Eqs. (13.46), (13.51), (13.55) and (13.56)). 

Before proceeding to a quadratic algebra approach to the problem in Section 3, it is 
worth making a few valuable observations. 

Mathematically speaking, the separable Schrodinger equation of our model admits four 
linearly independent solutions obtained by combining the two independent solutions of the 
second-order differential equation in x with those of the second-order differential equation 
in y. Among those four functions, only the combination ip^j(x,y), considered in (12.121) . 
satisfies all the boundary conditions and is normalizable onD. It is indeed clear that the 
alternative solution to the differential equation in x is not normalizable, while that to the 
differential equation in y, 

f sin[(/ + l)gy] for/ = 0,2,4,..., 

Xi{y) < (2-18) 
I cos[(Z + l)qy] for I — — 1, 1, 3, 5, . . . , 

violates the second condition in Eq. (12. 2p . Hence the three remaining combinations provide 
unphysical functions. 
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Some mathematical considerations might also lead to another choice than L and R^ k ' 
for the basic integrals of motion complementing H^ k \ First of all, instead of L, one might 
select the operator p y = —id y , which obviously satisfies the condition [H^ k \p y ] = 0. This 
would result in a linear and a quadratic (in the momenta) integrals of motion, generating 
a much simpler quadratic algebra than that to be considered in Section 3. It should be 
realized, however, that the eigenf unctions e imy (m e Z) of p y , being linear combinations 
of the physical and unphysical functions (12.141) and (12.181) . are useless from a physical 
viewpoint. We are therefore forced to consider the second-order operator L instead of p y . 

Furthermore, it is straightforward to see that another pair of first-order differential 
operators 

fj(k)j _ _ cos \y g X cos gy Q x _ gjnh g X s j n qyQ^_q ginh qx cos qy 

— qk csch qx cos qy , (2-19) 
f]^ = cosh qx cos qy d x + sinh qx sin qy d y + q sinh qx cos qy 

— qk csch qx cos qy, (2.20) 

intertwines with H^ k ' and H\(k), i.e., satisfies the relation 

and its Hermitian conjugate. Such operators correspond to the choice a = c = g = 0, 
b = d = 1 in Eq. (12.29). 

As a consequence of (12.211) . the operator 

R (k) = ^(fc)t^) 

= — cosh 2 qx cos 2 qyd\ — 2 sinh qx cosh qx sin qy cos qy d 2 xy — sinh 2 qx sin 2 qy d y 
+ q sinh qx cosh qx(l — 4 cos 2 qy)d x — q(l + 4 sinh 2 gx) sin qy cos g?/<9 y 
+ g 2 (sinh 2 qx — cos 2 qy — 3 sinh 2 cos 2 gy) — g 2 fc(l + csch 2 qx cos 2 qy) 
+ q 2 k 2 csch 2 gx cos 2 qy, (2.22) 

commutes with and is therefore another integral of motion. It can of course be ex- 
pressed in terms of H^ k \ L and R^ k \ as it can be checked that 

H {k) = L + R {k) + R {k) + 2q 2 k. (2.23) 
7 



However, we have now at our disposal three (dependent) integrals of motion L, and 
in addition to H^ k \ so that we may ask the following question: what is the best choice 
for the basic integrals of motion from a physical viewpoint? 

This problem is easily settled by noting that the zero modes of ff k \ 

uf\x,y) = (tanh qx) k (sech qx) s+1 (sin qy) s , (2.24) 

violate the second condition in Eq. (12.21) for any real value of s and therefore lead to 
unphysical functions. This contrasts with what happens for the zero modes uis\x,y) of 
rf- k \ given in (13.28), which are physical functions for s > and can therefore be used to 
build the functions ty^ No (x, y) considered in (12.16p . as it was shown in (13.32). We conclude 
that the physics of the model imposes the choice of L and R^ as basic integrals of motion. 

3 Quadratic algebra approach to the PDM model in 
a two-dimensional semi-infinite layer 

3.1 Quadratic associative algebra and its classical limit 

It has been shown [221 12H] that for any two-dimensional quantum superintegrable system 
with integrals of motion A, B, which are second-order differential operators, one can con- 
struct a quadratic associative algebra generated by A, B, and their commutator C. This 
operator is not independent of A, B, but since it is a third-order differential operator, it 
cannot be written as a polynomial function of them. The general form of the quadratic 
algebra commutation relations is 

[A,B] = C, (3.1) 
[A,C] = aA 2 + ^{A,B} + 5A + eB + (, (3.2) 
[B,C] = aA 2 - 7 5 2 -a{A,B} + dA-5B + z. (3.3) 

Here {A, B} = AB + BA, 

6 = 6(H)^S + 6 1 H, e = e(H)=e + e 1 H, ( = ((H) = Co + Cl# + C2H 2 , 
d = d(H) = d + d 1 H 1 z = z(H) = z + z x H + z 2 H 2 , (3.4) 

and a, 7, a, Si, e«, Q, d iy Z\ are some constants. Note that it is the Jacobi identity 
[A, [B, C]] = [B, [A, C]] that imposes some relations between coefficients in (I3.2p and (13.31) . 
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Such a quadratic algebra closes at level 6 [28] or, in other words, it has a Casimir 
operator which is a sixth-order differential operator [22J, 

K = C 2 + laA i -\a{A,A 1 B}-\ 1 {A 1 B 1 B}+{la 2 + d+la 1 )A 2 
+ (fary - 5) {A, B} + (| 7 2 - e) £? 2 + (forf + |ae + ±d 7 + 2z) A 
+ (-|ae + §7<J - 2C) 5 + | 7 z - §c< 
= ko + k.H + h^ + hH 3 , (3.5) 

where ^ are some constants and {A, B, C} = ABC + ACS + £MC + BCM + + A 
For our two-dimensional PDM model, described by the Hamiltonian defined in Eqs. {\2.3\i 
- ( 12. 51) . we shall take 

A = R, B = L, (3.6) 

where, for simplicity's sake, we dropped the superscript (k) because no confusion can arise 
outside the SUSYQM context. 

To determine their commutation relations, it is worth noting first that their building 
blocks, the first-order differential operators d y , rj* and r], generate another quadratic algebra 
together with the other set of intertwining operators fy , fj, given in f 1 2 . X 9 j) and (12.201) . Their 
commutation relations are indeed easily obtained as 

[d y ,v] = QV, [9 v ,fj] = -<H, [v,v} = qd y , (3.7) 
[77,77+] = 2q 2 k(l + e), [vrf] = 2q 2 k(l+?), [v, tf] = -qd y + 2q 2 k& (3.8) 

and their Hermitian conjugates. In (13.81) . we have defined 

£ = —(2qk)~ 1 (rj + rj') = csch qx sin qy, £ = — {2qk)~ 1 {fj + fj*) = csch qx cos qy. (3.9) 

Interestingly, d y , i] and fj (as well as d y , rj* and fj*) close an sl(2) subalgebra. 

From these results, it is now straightforward to show that the operator C in (13 .ip is 
given by 

C = q{d y ,r]*fj + fj*r)} (3.10) 
and that the coefficients in (13.21) and (13. 3p are 

a = 7 = 8g 2 , S = 8q 2 [q 2 (2k - 1) - H], e = 16q 4 (k - l)(k + 1), 

( = 8q\k- l)(2q 2 k-H), a = 0, d = lQq\ z = 8q\2q 2 k - if). (3.11) 
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On inserting the latter in (13. 5p . we obtain for the value of the Casimir operator 

K = -Aq A [2q 2 (7k - 6) - 3H](2q 2 k - H). (3.12) 

It is worth noting that since a = in (13.31) . we actually have here an example of quadratic 
Racah algebra QR(3) [19]. 

Before proceeding to a study of its finite-dimensional irreducible representations in Sec- 
tion 3.2, it is interesting to consider its classical limit. For such a purpose, since we have 
adopted units wherein h = 2m = 1, we have first to make a change of variables and of 
parameters restoring a dependence on h (but keeping 2m = 1 for simplicity's sake) before 
letting h go to zero. 

An appropriate transformation is 

X = hx, Y = hy, P x = -ihd x , Py = -ifrd Y , Q = f , K = hk. 

n 

(3.13) 

On performing it on the Hamiltonian given in Eqs. (12. 3p - (12. 5p . we obtain 

H = -h 2 (d x cosh 2 QXd x + d Y cosh 2 QXdy) - h 2 Q 2 cosh 2 QX + Q 2 K(K - h) csch 2 QX, 

(3.14) 



yielding the classical Hamiltonian 



H c = lim H = cosh 2 QX(P 2 + P 2 ) + Q 2 K 2 csch 2 QX. (3.15) 
h— >o 

A similar procedure applied to the intertwining operators leads to 
?7 r = lim 77 

= i cosh QX sin QYP X - i sinh QX cos QYPy - QK csch QX sin QY, (3.16) 
7/v = lim 7? 

= i cosh QX cos QYP X + i sinh QX sin QYP Y - QK csch QX cos QY, (3.17) 

together with r/* = lim^o 7 ?^ an d ffc = lim^o^^, while the operators quadratic in the 
momenta give rise to the functions 

L c = limL = P Y , (3.18) 

h— >o 

R c = UmR = cosh 2 QX sin 2 QYP 2 - 2 sinh QX cosh QX sin QY cos QYP X P Y 

h— >o 

+ sinh 2 Q A cos 2 QYPy + Q 2 A 2 csch 2 QX sin 2 QF, (3.19) 

i? c = lim R = cosh 2 QX cos 2 QFP|- + 2 sinh QX cosh QA sin QY cos QY P X P Y 
h— >o 

+ sinh 2 QA sin 2 QFPy + Q 2 A 2 csch 2 QA cos 2 QY, (3.20) 

10 



satisfying the relation 

H C = L C + R C + R c . (3.21) 

The quadratic associative algebra (13. ip - (13.51) is now changed into a quadratic 
Poisson algebra, whose defining relations can be determined either by taking the limit 
lim^ (i^) _1 [O, O'} = {O c , 0' c }p or by direct calculation of the Poisson brackets {O c , 0' c }p: 

{A C ,B C } P = C c , (3.22) 
{A C ,C C } P = a c A 2 c + 2 lc A c B c + 5 c A c + e c B c + ( c , (3.23) 
{B C ,C C } P = a c A\ - lc B 2 c - 2a c A c B c + d c A c - 5 C B C + z c . (3.24) 



Here 



C c = lim % = 2QP Y (r ] * c r ]c + fj* Vc ) (3.25) 
ft— »o in 



and 

« c = 7c = -8Q 2 , 5 C = 8Q 2 H C , e c = -16Q 4 K 2 , ( c = a c = d c = z c = 0. (3.26) 
Such a Poisson algebra has a vanishing Casimir: 

K c = lim K = 0. (3.27) 

h— >o 

3.2 Finite-dimensional irreducible representations of the 
quadratic associative algebra 

The quadratic algebra (13. ip - (13.51) can be realized in terms of (generalized) deformed 
oscillator operators Af, b\ b, satisfying the relations 



[J\f, 6 f ] = b\ [J\f, b] = -b, tfb = $(A0, btf = $(AT + 1), (3.28) 
where the structure function is a 'well-behaved' real function such that 

$(0) = 0, $(x) > for x > 0. (3.29) 
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m = 0,1,2, 



(3.30) 



This deformed oscillator algebra has a Fock-type representation, whose basis states |m), 
m = 0, 1, 2, ... [j] fulfil the relations 

Af\m) = m\m), 

tf\m) = y/&(m+ 1) \m + 1), 
6|0) = 0, 

b\m) = a/ $(m) | m — 1), 771 = 1, 2, ... . 
We shall be more specifically interested here in a subclass of deformed oscillator opera- 
tors, which have a (p + l)-dimensional Fock space, spanned by \p, m) = \m), m = 0, 1, . . . , 
p, due to the following property 

$(p+l)=0 (3.31) 
of the structure function, implying that 

(&t)P+i = h p+l = 0. (3.32) 

These are so-called (generalized) deformed parafermionic oscillator operators of order p |29j . 
The general form of their structure function is given by 



= x(p + 1 — x)(a + aix + a 2 x 2 + • ■ • + a p _ix p 1 ] 



(3.33) 



where ao, a%, . . . , a p _i may be any real constants such that the second condition in f!3.29j) 
is satisfied for x — 1, 2, . . . , p. 

A realization of the quadratic algebra (13.11) - f)3.5p in terms of deformed oscillator 
operators TV, tf, b reads 



A = A(J\f), (3.34) 
B = a{M) + 6 f p(A0 + p(Af)b, (3.35) 

where A(Af), cr(Af) and p(Af) are some functions of Af, which, in the 7 ^ case, are given 
by 



A(Af) 
a(AT) 
p\Kf) 



7 
2 



(Af + u) 



1 e 
4 ~ ~ 2 



a 
"I 



(A/" + M ) 



+ 



ae — 7<5 ae 2 — 2^5e + 47 2 ( 
2 7 2 



4 7 4 



(AT + M ) 



3 ■ 2 12 7 8 (AT + u )(jV + u + 1) [2(JV + u) + 1] : 



(3.36) 

T .(3-37) 
4 

(3.38) 



x We adopt here the unusual notation \m) in order to avoid confusion between the number of deformed 
bosons and the quantum number n introduced in (|2. 12[) . 
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with the structure function 



= -3072 7 6 ir[2(7V + M) - l] 2 

- 48 7 6 (a 2 e - a 7 <5 + a 7 e - d 7 2 ) [2 (TV + u) - 3] [2 (AT + u) - 1] 4 [2(A/" + u) + 1] 
+ 7 8 (3a 2 + 4a 7 )[2(A/" + w) - 3] 2 [2(A/" + it) - l] 4 [2(Af + u) + l] 2 

+ 768(ae 2 - 2 7 5e + 4 7 2 () 2 

+ 32 7 4 (3a 2 e 2 - 6a 7 5e + 2a 7 e 2 + 2 7 2 5 2 - 4d 7 2 e + 8 7 3 z + 4a 7 2 () 
x [2(M + u) - \\ 2 \\2{U + u) 2 - \2{U + u) - 1] 

- 256 7 2 (3aV - 9a 7 5e 2 + a 7 e 3 + 6 7 2 5 2 e - 3d 7 V + 2 7 4 <5 2 + 2c/ 7 4 e + I2^ez 



These functions depend upon two (so far undetermined) constants, u and the eigenvalue of 
the Casimir operator K (which we denote by the same symbol). 

Such a realization is convenient to determine the representations of the quadratic alge- 
bra in a basis wherein the generator A is diagonal together with K (or, equivalently, H) 
because the former is already diagonal with eigenvalues A(m). The (p + l)-dimensional 
representations, associated with (p+ l)-fold degenerate energy levels, correspond to the re- 
striction to deformed parafermionic operators of order p [22]. The first condition in f!3.29j) 
can then be used with Eq. (13.311) to compute u and K (or E) in terms of p and the Hamil- 
tonian parameters. A choice is then made between the various solutions that emerge from 
the calculations by imposing the second restriction in (13.291) for x — 1, 2, . . . , p. 

In the present case, for the set of parameters ( 13. Ill) , the complicated structure function 



(13.391) drastically simplifies to yield the factorized expression 

$(x) = 3-2 3 y°(2x + 2w + £;- l)(2x + 2w + A; - 2)(2x + 2w - k)(2x + 2u-k- 1) 
x (2x + 2u - \ + A) [2x + 2u - § + A) (2x + 2u - \ - A) 



4 7 5 z + 12a 7 2 eC - 12 7 3 5( + 4a 7 4 () [2(JV + u) - if. 



(3.39) 



x (2x + 2u 



3 
2 



A), 



(3.40) 



where 




(3.41) 



Furthermore, the eigenvalues of the operator A become 



A(m) = q 2 (2m + 2u - k)(2m + 2u + k). 



(3.42) 
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Since A = R is a positive-definite operator, only values of u such that A (to) > for m = 0, 
1, . . . , p should be retained. 

On taking this into account, the first condition in (I3.29|) can be satisfied by choosing 
either u = k/2 or u = (k + l)/2, yielding 

A(m) = Aq 2 m(m + k) (3.43) 

or 

A{m) = Aq 2 (m + §) (m + k + §) , (3.44) 

respectively. For u = k/2, Eq. ( 13. 3 1 f) . together with the second condition in (13.291) . can be 
fulfilled in two different ways corresponding to A = 2p + k + 1 ± | or 

E = q 2 (2p+\±\) (2p + 2k + \±\). (3.45) 

The resulting structure function reads 

$(x) = 3 • 2 3 y°xG» + 1 - x) (x - i) (p+l±i-x) (z + fc-i) (z + fc-l) 

x(x + p + A; + i±i)(x + p + A;-i±i). (3.46) 

Similarly, for tt = (A; + l)/2, we obtain 

E = q 2 (2p + | ± \) (2p + 2k + | ± |) (3.47) 

and 

$(x) = 3 • 2 38 q 20 x{p +l-x)(x + \)(p + l±\-x){x + k)(x + k-\) 

x {x+p + k + \±\){x + p + k + \±\) . (3.48) 

Our quadratic algebra approach has therefore provided us with a purely algebraic deriva- 
tion of the eigenvalues of H and R in a basis wherein they are simultaneously diagonal. It 
now remains to see to which eigenvalues we can associate physical wavefunctions, i.e., nor- 
malizable functions satisfying Eq. (12. 2\\ . This will imply a correspondence between \p, m) 
and the functions tyjy t N- u (x, y), defined in (12.161) . 

On comparing A(m) to the known (physical) eigenvalues r v of R, given in (I2.17p . we 
note that the first choice (13.431) for A(m) corresponds to even v = 2m (hence to even N), 
while the second choice (I3.44p is associated with odd v = 2m + 1 (hence with odd N). 
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Appropriate values of p leading to the level degeneracies (12.111) are therefore p = N/ 2 and 
p = (N — l)/2, respectively. With this identification, both Eqs. (I3.45P and (I3.47P yield the 
same result 

E = q 2 {N +l±\) (N + 2k + \±\) . (3.49) 

Comparison with (12. lOj) shows that only the upper sign choice in (I3.49p leads to physical 
wavefunctions ^n,n-u(x, y). 

Restricting ourselves to such a choice, we can now rewrite all the results obtained in 
this subsection in terms of iV and v instead of p and m. In particular, the two expressions 
(13.461) and (13.481) for the structure function can be recast in a single form <3>(m) — > 
where 

= 3-2™q 2Q u(u-l)(u+2k-l)(u+2k-2)(N +u+2k)(N +u+2k + l)(N -u+2)(N -u+3). 

(3.50) 

More importantly, our quadratic algebraic analysis provides us with an entirely new 
result, namely the matrix elements of the integral of motion L in the basis wherein H and 
R are simultaneously diagonal. On using indeed the correspondence \p, m) — > ^n^_ u , as 
well as the results in Eqs. ( 13301 . ([335]) . (13371) . fl33BD and (13301 . we obtain 

L^N,N-u = &v^N,N~v + T v ^ N ^ N - v+ 2 + T u+2 ^ N,N-is-2, (3.51) 

where we have reset a{m) — > a u , p{m) — > p u and defined r u = s v p v ^2\P&l- The explicit 
form of the coefficients on the right-hand side of (I3.5ip is given by 

+ [N 2 + (2k + 3)N + 2k 2 + 2k + + k - l)(y + k + 1) 

-k(k-l)(N + k + l)(N + k + 2)}, (3.52) 

x (N - v + 2)(N - v + 3)(N + v + 2k){N + v + 2k + 1). (3.53) 

Note that tv is determined up to some phase factor s v depending on the convention adopted 
for the relative phases of ^n,n-v and ^n,n-v+2- 
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For N = 4, for instance, v runs over 0, 2, 4, so that Eqs. (I3.5ip - (13.531) become 
q 2 \, x l2(k+l)(2k + 3)(2k + 9) 



L\I> 



4,2 



A; + 3 

3^4 



k + 2 



1.2 



(3.54) 



L\I> 



+ 



q 



1:1 



fc + 3 

■§2 

k + 1 

2 

«2 



2(/c + l)(2k + 3)(2A; + 9) 17fc 2 + 76k + 39 T 
-# 4 o + — — — # 



k + 2 



(k + l)(k + 3) 



4,2 



10(A; + 3)(2fc + l)(2fc + 7) 



k + 2 



i,-i 



10(fc + 3)(2£; + l)(2£; + 7) 



k + 2 



#4,2 + 5(fc + 3)^4,4 



(3.55) 
(3.56) 



As a check, these results can be compared with those derived from the action of L on the 
expansions of #4,2 and ^4 4 in terms of the first basis eigenf unctions ^ ,4 5 ^1,2 and -02,o 
(see, e.g., Eqs. (13.61) and (13.49) for ^0 and #4,4, respectively). This leads to the phase 
factors s 2 = s 4 = — 1. 

To conclude, it is worth mentioning that had we made the opposite choice in Eq. (13.61) . 
i.e., A = L and B = R, we would not have been able to use the deformed parafermionic 
realization (I3.34p . (I3.35P to determine the energy spectrum. The counterpart of the 
parafermionic vacuum state would indeed have been a function annihilated by L and there- 
fore constructed from the unphysical function X-iiv) °f Eq. (12.181) . 

4 Exactly solvable PDM models in three dimensions 



In the present section, we plan to show that the Hamiltonian ( 12.41) on the two-dimensional 
domain (12.11) can be easily extended to three dimensions in such a way that the domain keeps 
its essential characteristic of abrupt termination while the Hamiltonian remains exactly 
solvable. The latter will still be integrable with three independent integrals of motion 
H, L and M, but the superintegrability of the two-dimensional model will be lost. This 
generalization can be carried out in two different ways. 

4.1 Exactly solvable PDM model in a semi-infinite paral- 
lelepipedal channel 

In (12.41) . let us replace the operator d 2 by the two-dimensional Laplacian in cartesian co- 
ordinates d 2 + &l and assume that z varies in the same range as y. This leads us to the 
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Hamiltonian 



H = —d x cosh 2 qxd x — cosh 2 qx(d 2 + d 2 ) — q 2 cosh 2 qx + q 2 k(k — 1) csch 2 qx, (4.1) 

defined on the semi-infinite parallelepipedal domain 

D: 0<x<oo, ~ <y,z < (4.2) 

2q 2q 

with wavef unctions satisfying the conditions 

V>(0,j/,*) = 0, ,l>(x,±^,zj=0, V (x, y, ±|J = 0. (4.3) 

Such a Hamiltonian commutes with the operators 

L = -d 2 , M = -d 2 z . (4.4) 

Their simultaneous normalizable eigenf unctions if) n> i tm (x,y,z), fulfilling (14.31) . can be easily 
obtained along the lines detailed in Section 3.1 of I. They can be written as 

(x, y, z) = (f) n ,i,m(x)Xi(y)(m(z), (4-5) 

where Xl(v) is given in (12.141) . ( m (z) can be expressed in a similar way with m and z 
substituted for I and y, respectively, while 

(f>n,i,m( x ) = A/;, iim (tanhgx) fe (sechgx) 1+d P n v J (1 - 2tanh 2 qx), (4.6) 

with 



5 = ^(l + iy + (m + iy, (4.7) 

/ 2g (2n + k + l + 5)n\T(n + k + ± + 5) \ 1/2 
^ " ^ T(n + l + 5)T(n + k + ±) ) ' ^ 

The simultaneous eigenvalues of L, M and H are (/ + l) 2 q 2 , (m + l) 2 q 2 and 

^n,i, m = g 2 (2n + 1 + 5) (2n + 2A; + 5) , (4.9) 

where n, i, m — 0, 1, 2, . . . . The only remaining degeneracies are those connected with the 
(l,m) exchange, i.e., E ny i )m = E n)m> \ for I ^ m, as well as some 'accidental' degeneracies, 
such as -E n ,i,8 = E n £fi corresponding to 5 = \/85. 
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4.2 Exactly solvable PDM model in a semi- infinite cylindrical 
channel 

Alternatively, we may replace d 2 in (j2.4p by the two-dimensional Laplacian in polar coor- 
dinates d 2 + - p d p + -^d 2 and assume that p varies in a finite domain, on the boundary of 
which wavefunctions vanish. In this way, we get the Hamiltonian 

H = —d x cosh 2 qxd x — cosh 2 qx ^d 2 H — d p H — - cf^j — q 2 cosh 2 qx + q 2 k{k — 1) csch 2 qx, 

(4.10) 

defined on the semi-infinite cylindrical domain 

D : 0<x < oo, < p< R, < ^ < 2vr, (4.11) 
with wavefunctions such that 

</>(0,p,<p) =0, i/>(x,R,<p) = 0, *P{x,p,2<K)=iP(x,p,0). (4.12) 
The two operators commuting with H are now 

L = -f<9 2 + ~^ + i<9 2 j, M = -id v . (4.13) 
The simultaneous normalizable eigenfunctions of H, L and M can be written as 

Here 

Uip) = (4.15) 
y lit 

corresponds to the eigenvalues m — 0, ±1, ±2, ... of M. Furthermore, 

X\m\,s(P) =Af\ m \, a J\m\(K\m\,aP), K \m\,s = (4.16) 

where J| m |(z) is a Bessel function, the symbol j\ m \ s , s = 1, 2, . . . , conventionally denotes [3T] 
its real, positive zeros, and [32] 

A/] m | iS = v^2 [RJ H+1 {j lmls )] _1 , (4.17) 

provides normalized solutions to the eigenvalue equation 

(p)Cm(<P), (4-18) 
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which satisfy the second condition in (l4.12p Fl Finally, 4> n ,\m\,s{.x) and the energy eigenvalues 
E n ,\m\,s are still given by the right-hand sides of Eqs. (14. 6[) and (14.91) . but with S now defined 
by 

s = K\m\,s = j\m\,B ^ , 41g s 

q qR 

This time the only level degeneracy left is that connected with the sign of m. 

5 Conclusion 

In this paper, we have revisited the exactly solvable PDM model in a two-dimensional 
semi-infinite layer introduced in I. Here we have taken advantage of its superintegrability 
with two integrals of motion L and R that are quadratic in the momenta to propose a third 
method of solution in the line of some recent analyses of such problems. 

We have first determined the explicit form of the quadratic associative algebra generated 
by L, R and their commutator. We have shown that it is a quadratic Racah algebra QR(3) 
and that its Casimir operator K is a second-degree polynomial in H. We have also obtained 
the quadratic Poisson algebra arising in the classical limit. 

We have then studied the finite-dimensional irreducible representations of our algebra 
in a basis wherein K (or H) and R are diagonal. For such a purpose, we have used a simple 
procedure, proposed in [22], consisting in mapping this basis onto deformed parafermionic 
oscillator states of order p. Among the results so obtained for the energy spectrum, we have 
selected those with which physical wavefunctions can be associated. This has illustrated 
once again the well-known fact that in quantum mechanics the physics is determined not 
only by algebraic properties of operators, but also by the boundary conditions imposed on 
wavefunctions. Our analysis has provided us with an interesting new result, not obtainable 
in general form in the SUSYQM approach of I, namely the matrix elements of L in the 
basis wherein H and R are simultaneously diagonal. 

In the last part of our paper, we have extended our two-dimensional model to three 
dimensions in two different ways by considering either a semi-infinite parallelepipedal chan- 
nel or a semi-infinite cylindrical one. Both resulting models remain integrable and exactly 
solvable, but the superintegrability of the two-dimensional model is lost. From a physical 



2 For the solution of a similar problem see [33] , 
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viewpoint, they illustrate the generation of bound states in a quantum channel when the 
uniformity is broken by an abrupt termination. 

As a final point, it is worth observing that the procedure used here to construct the 
irreducible representations of QR(3) is not the only one available. In particular, the ladder- 
operator method employed in [T9j [20] would allow us to express the transformation matrix 
elements between the bases ip^i an d *nn (denoted by , in I) in terms of Racah- 

Wilson polynomials. 
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